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Abstract — The Alpha Magnetic Spectrometer (AMS), whose 
final version AMS-02 is to be installed on the International 
Space Station (ISS) for at least 3 years, is a detector designed to 
measure charged cosmic ray spectra with energies up to the TeV 
region and with high energy photon detection capability up to 
a few hundred GeV, using state-of-the art particle identification 
techniques. 

Among several detector subsystems, AMS includes a proximity 
focusing RICH enabling precise measurements of particle electric 
charge and velocity. The combination of both these measurements 
together with the particle rigidity measured on the silicon tracker 
endows a reliable measurement of the particle mass. 

The main topics of the AMS-02 physics program include 
detailed measurements of the nuclear component of the cosmic- 
ray spectrum and the search for indirect signatures of dark 
matter. Mass separation of singly charged particles, and in 
particular the separation of deuterons and antideuterons from 
massive backgrounds of protons and antiprotons respectively, is 
essential in this context. Detailed Monte Carlo simulations of 
AMS-02 have been used to evaluate the detector's performance 
for mass separation at different energies. The obtained results 
and physics prospects are presented. 

I. The AMS-02 experiment 

The Alpha Magnetic Spectrometer (AMS)[1], whose final 
version AMS-02 is to be installed on the International Space 
Station (ISS) for at least 3 years, is a detector designed to study 
the cosmic ray flux by direct detection of particles above the 
Earth's atmosphere using state-of-the-art particle identification 
techniques. AMS-02 is equipped with a superconducting mag- 
net cooled by superfluid helium. The spectrometer is com- 
posed of several subdetectors: a Transition Radiation Detector 
(TRD), a Time-of-Flight (TOF) detector, a Silicon Tracker, 
Anticoincidence Counters (ACC), a Ring Imaging Cerenkov 
(RICH) detector and an Electromagnetic Calorimeter (ECAL). 
Fig. [T] shows a schematic view of the full AMS-02 detector. A 
preliminary version of the detector, AMS-01, was successfully 
flown aboard the US space shuttle Discovery in June 1998 [2]. 

The main goals of the AMS-02 experiment are: 

• A precise measurement of charged cosmic ray spectra 
in the rigidity region between ~ 0.5 GV and ~ 2 TV, 
and the detection of photons with energies up to a few 
hundred GeV; 

• A search for heavy antinuclei (Z > 2), which if discov- 
ered would signal the existence of cosmological antimat- 
ter; 
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Fig. 1 . Exploded view of the AMS-02 detector. 

• A search for dark matter constituents by examining 
possible signatures of their presence in the cosmic ray 
spectrum. 

The long exposure time and large acceptance (0.5 m 2 sr) of 
AMS-02 will enable it to collect an unprecedented statistics 
of more than 10 10 nuclei. 

II. The AMS RICH detector 

One of the subdetectors in AMS-02 is a proximity focusing 
Ring Imaging Cerenkov (RICH) detector. It is composed of 



a dual radiator with silica aerogel (n = 1.050) and sodium 
fluoride (n = 1.334), a high reflectivity lateral conical mirror 
and a detection matrix with 680 photomultipliers coupled to 
light guides. 
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Fig. 2. A simulated proton event as seen in the AMS-02 display. 
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Fig. 3. The same event of Fig. [2] as seen in the RICH display developed at 
LIP. 



the conical reflector was included in the beam test setup[5]. 
The algorithms for velocity and charge reconstruction were 
successfully applied to data from these prototype tests [6]. 



The RICH detector will provide a very accurate velocity 
measurement (in aerogel, Af3/(3 ~ 10~ 3 and 10~ 4 for Z = 1 
and Z = 10 — 20, respectively) and charge identification of 
nuclei up to iron (Z = 26). 

RICH data, combined with information on particle rigidity 
from the AMS Silicon Tracker, enable the reconstruction of 
particle mass. A typical RICH event is shown in Figs. [2] and 
[3] where the latter gives a detailed view of the readout matrix. 
The accuracy of the RICH velocity measurement is essential 
due to the growth of relative errors when v — ► c: 

Am Ap , A/3 

= — ©7 2 -f 

m p p 

The assembly of the AMS RICH detector is currently 
underway at CIEMAT in Madrid. The integration of the RICH 
and the other subdetectors of AMS-02 will take place at CERN 
in 2007. 

The analysis of RICH data involves the identification of the 
Cerenkov ring in a hit pattern which usually includes several 
scattered noise hits and an eventual strong spot in the region 
where the charged particle crosses the detection plane. Two 
independent algorithms for velocity and charge reconstruction 
have been developed in the AMS collaboration for the analysis 
of RICH events: a geometrical method based on a hit-by-hit 
reconstruction^], and a method using all the hits with the 
maximization of a likelihood function[4]). 

A prototype of the RICH detector, consisting of 96 photo- 
multiplier units, was tested both with cosmic ray particles and 
with beam ions at the CERN SPS in 2002 and 2003. A piece of 



III. Components of the cosmic ray spectrum 

Protons are the most abundant component (~ 90%) of 
charged cosmic rays reaching the Earth's vicinity. The remain- 
ing fraction is essentially made of atomic nuclei (mostly 4 He) 
and single-charged particles such as e~, e + and p. 

Precise measurements of the smaller components in the 
cosmic ray spectrum are essential in the context of the study 
of cosmic -ray production and acceleration. The required pre- 
cision can only be attained through very effective charge and 
mass discrimination methods since the abundances of different 
components differ by several orders of magnitude. 

Ratios such as D/p, 3 He/ 4 He and B/C give information on 
the interstellar medium since all compare the abundances of 
secondary and primary species. The beryllium isotope ratio 
10 Be/ 9 Be is a probe for galactic confinement times since both 
isotopes are secondary but one of them, 10 Be, is unstable 
(t 1/2 = 1.5 x 10 6 yr). 

The lightest neutralino (xj), predicted by supersymmetric 
models, is a strong dark matter candidate. If it exists and 
accounts for the known dark matter density (CIcdm — — 
57b ~ 0.2) or at least for a significant part of it, neutralino 
annihilation (xiXi ~^ •••) must take place and contribute to 
the observed cosmic ray composition, with the more visible 
effects occurring in the spectra of antiparticles like e + , p 
and especially D[7]. Fig. |4] shows a comparison between the 
expected D fluxes from secondary production and from dark 
matter annihilation. 
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Fig. 4. Comparison between expected antideuteron flux from secondary 
processes (dashed line) and the flux from the annihilation of a 60 GeV dark 
matter particle (solid lines: solar minimum; dotted lines: interstellar flux). 
Fluxes from dark matter annihilation are shown for three sets of propagation 
parameters, (from Ref. [7]) 

IV. Particle identification 

To evaluate the capabilities of AMS-02 for mass separation 
of deuterons and antideuterons from other particles with the 
same charge, studies have been performed using the case of 
deuteron vs. proton separation. In the past, studies on the 
separation of helium (Z = 2) and beryllium (Z = 4) isotopes 
have also been performed using a standalone simulation of the 
RICH detector[8]. In the present case the large difference be- 
tween proton and deuteron abundances (D/p ~ 1%) increases 
the importance of a very effective mass separation to isolate 
the deuteron signal from a large background of proton events. 

In the study of D/p separation a full-scale simulation of the 
AMS detector was used. Particles were simulated as coming 
from the top plane of a cube, corresponding to an acceptance 
of 47.78 m 2 sr. Three data samples were chosen. TableUshows 
the momentum ranges and number of events simulated in each 
sample. 



TABLE I 

Samples used in the D/p separation studies 



Sample 


Momentum range 


No. events 


p (low momentum) 


0.5 - 10 GeV/c 


3.1 x 10 B 


p (high momentum) 


10 - 200 GeV/c 


1.3 x 10 s 


D 


0.5 - 20 GeV/c 


5.6 x 10' 



For each sample, g^pi = constant. Variable weights were 
assigned to events in order to compensate for the statistics in 
each sample and to reproduce a realistic spectrum (Fig. |5}: 

> The simulated proton spectrum followed 

dN/dE oc E~ 2 - 7 ; 



Kinetic energy (GeV/nucleon) 

Fig. 5. Simulated proton and deuteron spectra used in this work. 

• The simulated deuteron spectrum was calculated combin- 
ing the proton spectrum above with D/p ratios taken from 
Ref. [9]. 

In each event a set of preliminary data selection cuts using 
readings from different subdetectors of AMS-02 was applied 
to reduce the fraction of events with a bad reconstruction. Only 
downgoing events (f3 > 0) were accepted. In addition, events 
were accepted if the following conditions were satisfied: 

• Only one particle was detected in the event; 

• A particle track was reconstructed by the Silicon Tracker; 

• No clusters were found in the Anti-Coincidence Counters; 

• Clusters from at least 3 TOF planes (out of 4) were used 
for event reconstruction; 

• At most one additional cluster was allowed in the TOF; 

• At least 6 Tracker layers (out of 8) were used in the track 
reconstruction; 

• Compatibility was required for the rigidity measure- 
ments obtained from two different algorithms, with 

AR/R < 3%; 

• Compatibility was also required for the rigidity measure- 
ments obtained from each half of the Tracker (upper and 
lower), with AR/R < 50%; 

• The particle's impact point on the RICH radiator was less 
than 58 cm from the centre (i. e. more than 2 cm from 
the mirror); 

• At most one track was present in the TRD; 

• The TOF and Tracker charge reconstructions were com- 
patible. 

Among the events that triggered the detector, a fraction 
corresponding to ~ 15-20% of proton events and ~ 10-15% 
of deuteron events in the relevant region of kinetic energy (few 
GeV/nucleon) passed this set of preliminary cuts, correspond- 
ing to an acceptance of ~ 0.3 m 2 sr for protons and ~ 0.2 
m 2 sr for deuterons. 

The reconstruction of particle masses was then performed 
for events having a signal in the RICH detector. The ex- 
tremely accurate velocity measurement provided by the RICH 
{A/3/ (3 ~ 10~ 3 in the case of protons and deuterons) is crucial 
to reduce the background level. A series of event selection cuts 
were introduced, based on data provided by the RICH and the 



results of the two reconstruction algorithms: 
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Fig. 6. Examples of inverse mass distribution in aerogel events for two 
energy regions. 
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Fig. 7. Example of inverse mass distribution in NaF events. 

• A Cerenkov ring was reconstructed using each method, 
and at least 3 hits were used in both cases; 

> The total ring signal was not higher than 10 photoelec- 
trons in NaF events, and not higher than 15 photoelec- 
trons in aerogel events; 

• A Kolmogorov test to the uniformity of the hits azimuthal 
distribution in the ring gave a result of at least 0.2 in the 
case of NaF events, and 0.03 in the case of aerogel events; 
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Fig. 8. Expected sensitivity of AMS for D/p ratio with one day of data. 

• Compatibility was required for the velocity measurements 
from the TOF and RICH detectors, with A/3//3 < 10%; 

• Compatibility was also required for the velocity mea- 
surements obtained from the two RICH reconstruc- 
tion methods, with A/3//3 < 0.3% for NaF events, and 
A/3//3 < 0.1% for aerogel events; 

• The reconstructed, rounded electric charge obtained from 
the geometrical method was 1 or 2; 

• The reconstructed, non-rounded electric charge obtained 
from the likelihood method was between 0.5 and 1.5 in 
NaF events, and between 0.6 and 1.4 in aerogel events; 

• The ring acceptance (visible fraction), as estimated by the 
likelihood method, was at least 20% in NaF events, and 
at least 40% in aerogel events; 

• The number of noisy hits not associated to the crossing 
of the charged particle (i. e., hits that were far from the 
reconstructed ring and far from the estimated crossing 
point of the charged particle in the detection matrix) was 
not higher than 2 in NaF events, and not higher than 3 
in aerogel events. 

V. Analysis results 

Results show that mass separation of particles with Z = 1 
is feasible even if one species is orders of magnitude more 
abundant than the other. D/p separation is possible up to 
Ekin ~ 8 GeV/nucleon. Some examples of the mass distri- 
butions obtained are shown in Figs. [6] and [7] Solid lines show 
the mass distributions before the RICH cuts were taken into 
consideration. 

Fig. [8] shows the expected sensitivity of AMS for the D/p 
ratio after one day of data taking. Results show that a single 
day of AMS-02 statistics will be sufficient to improve on the 
existing data for this ratio. 



In the optimal region immediately above the aerogel ra- 
diation threshold (Ekin < 5 GeV/nucleon) rejection factors 
higher than 10 4 region were attained (Fig.|9]l. The best relative 
mass resolutions for protons (Fig. ITOb and deuterons are 
~ 2% for both radiators in the regions above their respective 
thresholds. 

After all cuts, an acceptance of ~ 0.07 m 2 -sr was ob- 
tained for protons, and ~ 0.05 m 2 -sr for deuterons at 
Ekin > 3 GeV/nucleon (Fig. fTTT i. The increase by a factor 
~ 10 in the acceptance above the aerogel threshold reflects the 
relative dimensions of the two radiators in the RICH detector. 

The main background in the deuteron case comes from non- 
gaussian tails of proton events with a bad velocity reconstruc- 
tion. Errors in rigidity reconstruction (AR/R ~ 2% in the 
GeV region) are not critical for this case. 

The specific set of cuts shown here corresponds to an ex- 
ample of a selection procedure. Other variations are possible. 
In particular, rejection factors may be improved by applying 
stricter cuts, at the expense of a further acceptance reduction. 
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Fig. 9. Rejection factor for D/p separation in aerogel events. 
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Fig. 10. Relative mass resolution for protons: NaF events (open dots) and 
aerogel events (filled dots). 
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Fig. 11. Acceptance for protons (top) and deuterons (bottom) at different 
stages of event analysis. The lower line in each plot corresponds to the final 
acceptance. 



10 10 events will be collected during its operation. Detailed 
simulations have been performed to evaluate the detector's 
particle identification capabilities, in particular those of the 
RICH. Simulation results show that the separation of light 
isotopes is feasible. Using a set of simple cuts based on 
event data, relative mass resolutions of ~ 2 % and rejection 
factors higher than 10 4 have been attained in D/p separation 
at energies of a few GeV/nucleon. The separation procedure 
presented here might be crucial for the identification of an 
antideuteron flux resulting from neutralino annihilation. 
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VI. Conclusions 



AMS-02 will provide a major improvement on the current 
knowledge of cosmic rays. A total statistics of more than 



